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1.  Introduction 


The  prediction  of  aerodynamic  coefficients  for  projectile  configurations  is  essential  in  the 
assessment  of  the  performance  of  new  designs.  Accurate  determination  of  aerodynamics  is 
critical  to  the  low-cost  development  of  new  advanced  guided  projectiles,  rockets,  missiles,  and 
smart  munitions.  Fins,  canards,  and  jets  can  be  used  to  provide  control  for  maneuvering 
projectiles  and  missiles.  The  flow  fields  associated  with  these  control  mechanisms  for  the  Army 
weapons  are  complex,  involving  three-dimensional  (3-D)  shock-boundary  layer  interactions,  jet 
interaction  with  the  free  stream  flow,  and  highly  viscous  dominated  separated  flow  regions 
[1, 2, 3].  The  jet  interference  can  extend  over  significant  portions  of  the  projectile  and  must  be 
modeled  correctly.  For  missiles,  jet  thrusters  have  been  studied  over  a  number  of  years  to 
provide  high-speed  aerodynamic  control.  These  thrusters  interact  with  the  surrounding  flow 
field,  and  the  resulting  jet  interaction  flow  field  again  is  complex.  Recently,  several  studies  have 
shown  that  tiny  micro-electro-mechanical  system  or  synthetic  unsteady  jets  can  significantly 
alter  the  flow  field  and  pressure  distributions  for  airfoils  and  cylinders  [4, 5, 6],  These  synthetic 
jets  are  active  control  devices  with  zero  net  mass  flux  and  are  intended  to  produce  the  desired 
control  of  the  flow  field  through  momentum  effects.  Many  parameters  such  as  jet  location,  jet 
velocity,  and  actuator  frequency  can  affect  the  flow  control  phenomenon. 

Smith  and  Glezer  [4]  have  conducted  an  excellent  study  of  the  flow  control  by  synthetic  jets  to 
provide  increased  fundamental  understanding  of  the  flow  physics.  Amitay  et  al.  [5] 
experimentally  investigated  flow  separation  control  on  a  cylinder  using  synthetic  jet  actuators. 
Their  work  showed  that  the  interaction  of  the  synthetic  jet  with  the  free  stream  flow  resulted  in  a 
virtual  modification  of  the  body  shape  and  significantly  increased  the  lift  force  as  a  result  of  the 
flow  reattachment.  Aerodynamic  flow  control  over  an  unconventional  airfoil  has  also  been 
demonstrated  by  Amitay  et  al.  [6]  to  enhance  post-stall  performance  with  actuators  operating  at 
frequencies  higher  than  the  characteristic  frequency  of  the  airfoil.  The  synthetic  jets  are  also 
being  investigated  for  possible  applications  to  improve  heat  transfer  and  drag  reduction  and  to 
enhance  mixing  [7]  in  combustors,  etc.  The  present  analysis  involves  these  synthetic  jets  for 
projectile  aerodynamic  control.  The  emphasis  in  the  present  research  is  to  provide  insight  into 
the  interaction  of  these  unsteady  jets  with  the  free  stream  flow  and  to  determine  the  feasibility  of 
these  jets  for  aerodynamic  control  of  a  subsonic  projectile. 

Computational  and  experimental  data  for  these  jet  interactions  are  very  limited.  Simple  theories 
cannot  predict  the  complex  flow  fields  associated  with  the  jet  interaction,  and  experimental  tests 
are  very  expensive.  To  help  reduce  experimental  costs,  computational  fluid  dynamics  (CFD)  are 
being  used  to  predict  these  complex  flows  and  to  provide  detailed  pressure  and  force  and 
moment  data.  There  have  been  several  recent  numerical  studies  [8,9]  of  flow  separation  control 
with  synthetic  jet  actuators.  He  and  Krai  [8]  have  used  Reynolds-averaged  Navier-Stokes 
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(RANS)  to  study  the  effect  of  jet  location  and  jet-forcing  frequency  on  the  lift  and  drag  forces  on 
an  airfoil.  The  jet  actuator  was  shown  to  increase  the  time-averaged  lift  but  also  increased  the 
amplitude  of  oscillation.  Lee  and  Goldstein  [9]  have  used  the  direct  numerical  simulation  (DNS) 
to  model  a  two-dimensional  (2-D)  synthetic  jet.  Although  this  simulation  produced  very  good 
results,  the  use  of  DNS  in  practical  3-D  flows  of  interest  is  prohibitive  because  of  its  computing 
resources  requirement.  Even  large  eddy  simulations  (LES)  [10],  in  which  large  eddies  are 
computed  directly  and  the  small  scales  are  modeled,  require  a  large  computational  cost  compared 
to  RANS  simulations.  While  the  RANS  method  works  well  for  steady  flows,  the  accuracy  of 
this  method  for  unsteady  flows  may  be  less  than  desired.  Since  the  large  energy-containing 
eddies  are  computed  by  the  LES  method,  this  technique  is  more  capable  of  handling  unsteady 
shear  layers  and  wakes,  etc.  Recently,  therefore,  a  hybrid  approach  [11, 12]  that  combines 
RANS  and  LES  has  been  developed  to  solve  practical  problems  of  interest  involving  unsteady 
flows  at  a  reasonable  computational  cost.  Both  RANS  and  hybrid  RANS-LES  models  have  been 
used  in  the  present  study. 

The  advanced  CFD  capability  used  here  solves  the  Navier-Stokes  equations  [13  through  16]  and 
incorporates  unsteady  boundary  conditions  for  the  simulation  of  the  synthetic  jets.  Also,  a 
hybrid  RANS-LES  turbulence  model  [12]  was  used  for  the  accurate  numerical  prediction  of 
unsteady  jet  flows.  Numerical  flow  field  computations  have  been  made  for  steady  and  unsteady 
jets  for  a  projectile  configuration  at  a  low  subsonic  speed. 

The  following  sections  describe  the  numerical  procedure  and  the  computed  results  obtained. 
Results  obtained  for  a  projectile  configuration  are  presented  at  Mach  0.1 1  and  several  angles  of 
attack,  0°  to  4°.  Computed  data  have  been  compared  with  experimental  data  provided  by  the 
Georgia  Tech  Research  Institute  (GTRI)  [IT]1. 


2.  Solution  Technique 


The  complete  set  of  3-D  time-dependent  Navier-Stokes  equations  is  solved  in  a  time-accurate 
manner  for  simulations  of  unsteady  jets.  A  commercially  available  code,  CFD++  [14, 15, 16],  is 
used  for  the  time-accurate  unsteady  CFD  simulations.  The  basic  numerical  framework  in  the 
code  contains  unified  grid,  unified  physics,  and  unified  computing  features.  The  user  is  referred 
to  these  references  for  details  of  the  basic  numerical  framework.  Here,  only  a  brief  synopsis  of 
this  framework  and  methodology  is  given. 

The  3-D,  time-dependent  RANS  equations  are  solved  by  the  following  finite  volume  method: 


l 


Private  communication  with  J.  McMichael,  Georgia  Tech  Research  Institute. 
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(1) 


J-JWrfF+|[F-G]-£M  =  jHrfF 

ot  y  j, 

in  which  W  is  the  vector  of  conservative  variables,  F  and  G  are  the  inviscid  and  viscous  flux 
vectors,  respectively,  H  is  the  vector  of  source  terms,  V  is  the  cell  volume,  and  A  is  the  surface 
area  of  the  cell  face. 

The  numerical  framework  of  CFD++  is  based  on  the  following  general  elements: 

1 .  unsteady  compressible  and  incompressible  Navier-Stokes  equations  with  turbulence 
modeling  (unified  physics); 

2.  unification  of  Cartesian,  structured  curvilinear,  and  unstructured  grids,  including 
hybrids  (unified  grid); 

3.  unification  of  treatment  of  various  cell  shapes  including  3-D  hexahedral,  tetrahedral, 
and  triangular  prism  cells,  2-D  quadrilateral  and  triangular  cells,  and  one-dimensional  linear 
elements  (unified  grid); 

4.  treatment  of  multi-block  patched  aligned  (nodally  connected),  patched  nonaligned, 
and  overset  grids  (unified  grid); 

5.  total  variation  diminishing  discretization  based  on  a  new  multi-dimensional 
interpolation  framework; 

6.  Riemann  solvers  providing  proper  signal  propagation  physics,  including  versions  for 
preconditioned  forms  of  the  governing  equations; 

7 .  consistent  and  accurate  discretization  of  viscous  terms  via  the  same  multidimensional 
polynomial  framework; 

8.  point- wise  turbulence  models  not  requiring  knowledge  of  distance  to  walls; 

9.  versatile  boundary  condition  implementation,  including  a  rich  variety  of  integrated 
boundary  condition  types  for  the  various  sets  of  equations;  and 

10.  implementation  on  massively  parallel  computers  based  on  the  distributed  memory 
message-passing  model  via  native  message-passing  libraries  or  message-passing  interface, 
parallel  virtual  machine,  etc.  (unified  computing). 

The  code  has  brought  together  several  ideas  about  convergence  acceleration  to  yield  a  fast 
methodology  for  all  flow  regimes.  The  approach  can  be  labeled  as  a  “preconditioned  implicit 
relaxation”  scheme.  It  combines  three  basic  ideas:  (1)  implicit  local  time  stepping, 

(2)  relaxation,  and  (3)  preconditioning.  Preconditioning  the  equations  ideally  equalizes  the 
eigenvalues  of  the  inviscid  flux  Jacobians  and  removes  the  stiffness  arising  from  large 
discrepancies  between  the  flow  and  sound  velocities  at  low  speeds.  The  use  of  an  implicit 
scheme  circumvents  the  stringent  stability  limits  suffered  by  their  explicit  counterparts,  and 
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successive  relaxation  allows  cells  to  be  revised  as  information  becomes  available  and  thus  aids 
convergence.  These  features  of  the  code  have  been  extremely  useful  in  the  present  numerical 
simulations  at  very  low  subsonic  speeds. 

The  CFD++  code  also  uses  a  multi-dimensional  interpolation  that  more  accurately  represents 
local  behavior  of  flow-dependent  variables.  Second  order  discretization  was  used  for  the  flow 
variables  and  the  turbulent  viscosity  equation.  The  turbulence  closure  is  based  on  topology- 
parameter-free  formulations.  Two-equation  and  higher  order  hybrid  RANS-LES  turbulence 
models  were  used  for  the  computation  of  turbulent  flows.  These  models  are  ideally  suited  to 
unstructured  bookkeeping  and  massively  parallel  processing  because  of  their  independence  from 
constraints  related  to  the  placement  of  boundaries  and/or  zonal  interfaces. 

For  computations  of  unsteady  jet  interaction  flow  fields  that  are  of  interest  here,  dual  time 
stepping  as  described  in  Section  2.1  was  used  to  for  time  accuracy.  In  addition,  special  jet 
boundary  conditions  were  developed  and  used  for  numerical  modeling  of  synthetic  jets. 

2.1  Dual  Time  Stepping 

CFD++  is  often  used  in  the  “dual  time-stepping  mode”  to  perform  transient  flow  simulations. 

The  term  “dual  time  step”  implies  the  use  of  two  time  steps.  The  first  is  an  “outer”  or  global  (and 
physical)  time  step  that  corresponds  to  the  time  discretization  of  the  physical  time  variation  term. 
This  time  step  can  be  chosen  directly  by  the  user  and  is  typically  set  to  a  value  to  represent  1/100 
of  the  period  of  oscillation  expected  or  forced  in  the  transient  flow.  This  time  step  is  applied  to 
every  cell  (not  separately  varying). 

An  artificial  or  “inner”  or  “local”  time  variation  term  is  added  to  the  basic  physical  equations. 
This  time  step  and  corresponding  “inner  iteration”  strategy  is  chosen  to  help  satisfy  the  physical 
transient  equations  to  the  desired  degree.  If  the  inner  iterations  converge,  then  the  outer  physical 
transient  equations  (their  discretization)  are  satisfied  exactly — otherwise,  approximately.  For  the 
inner  iterations,  the  time  step  is  allowed  to  vary  spatially.  Also,  relaxation  with  multi-grid 
(algebraic)  acceleration  is  employed  to  reduce  the  residues  of  the  physical  transient  equations.  It 
is  found  that  an  order  of  magnitude  reduction  in  the  residues  is  usually  sufficient  to  produce  a 
good  transient  iteration.  This  may  require  a  few  internal  iterations  to  achieve  (between  3  and  10, 
depending  on  the  magnitude  of  the  outer  time  step),  the  nature  of  the  problem,  the  nature  of  the 
boundary  conditions,  and  the  consistency  of  the  mesh  with  respect  to  the  physics  at  hand. 

2.2  Unsteady  Jet  Boundary  Condition 

CFD++  has  a  large  collection  of  boundary  conditions  (BC).  Each  BC  is  encoded  as  a  basic  form, 
along  with  a  collection  of  modifiers.  One  particular  BC  used  for  the  simulations  presented  here 
is  a  “oscillating  jet”  BC.  In  its  basic  form,  it  is  a  steady  inflow/outflow  BC  wherein  the  user 
supplies  the  velocity  normal  to  the  boundary  along  with  static  temperature  and  any  turbulence 
quantities.  When  the  velocity  provided  is  negative,  it  is  considered  to  be  an  inflow  and  when  it  is 
positive,  it  is  treated  as  an  outflow.  In  the  case  of  inflow,  the  static  temperature  and  turbulence 
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quantities  are  used,  along  with  the  inflow  velocity.  In  the  case  of  outflow,  only  the  velocity  is 
used.  At  inflow,  the  tangential  component  of  velocity  is  set  to  zero  and  at  outflow,  the  tangential 
component  is  extrapolated  from  the  interior.  At  outflow,  all  primitive  variables  except  normal 
velocity  are  extrapolated  from  the  interior.  At  inflow,  the  static  pressure  is  taken  from  the 
interior. 

The  first  modifier  available  for  this  BC  allows  the  velocity  to  oscillate.  The  base  velocity  is 
multiplied  by  an  amplitude  that  varies  as  sin  (2tift)  in  which  f  is  the  frequency  of  the  oscillation. 
Thus,  the  oscillating  velocity  can  cycle  from  being  positive  to  being  negative  and  back  within 
each  period  (or  from  being  negative  to  positive  and  back,  based  on  the  sign  of  the  input  for  the 
basic  BC  formulation). 

A  second  modifier  permits  the  steady  or  oscillating  inflow/outflow  to  be  on  over-certain  time 
intervals  and  off  during  other  intervals.  During  “on”  periods,  the  basic  or  the  basic  multiplied  by 
the  oscillating  amplitude  multiplier  (first  modifier)  is  used.  During  “off’  periods,  the  amplitude 
is  set  to  zero.  Therefore,  when  the  jet  velocity  is  zero,  the  boundary  condition  behaves  like  a  no¬ 
slip  condition. 

2.3  Hybrid  RANS-LES  Model 

Currently,  the  two  most  popular  forms  of  turbulence  closure,  namely,  ensemble-averaged  models 
(typically  based  on  the  RANS  equations)  and  LES  with  a  sub-grid-scale  model,  face  a  number  of 
unresolved  difficulties.  Specifically,  existing  LES  models  have  met  with  problems  related  to  the 
accurate  resolution  of  the  near- wall  turbulent  stresses.  In  the  near-wall  region,  the  foundations  of 
large  eddy  simulation  are  less  secure,  since  the  sizes  of  the  (anisotropic)  near- wall  eddies 
approach  those  of  the  Kolmogorov  scale,  requiring  a  mesh  resolution  approaching  that  of  a  direct 
numerical  simulation.  On  the  other  hand,  existing  ensemble-averaged  turbulence  models  are 
limited  by  their  empirical  calibration.  Their  representation  of  small  scale  flow  physics  cannot  be 
improved  by  refining  the  mesh,  and  over  short  time  scales,  they  tend  to  be  overly  dissipative  with 
respect  to  perturbations  around  the  mean,  often  suppressing  unsteady  motion  altogether. 

While  LES  is  an  increasingly  powerful  tool  for  unsteady  turbulent  flow  prediction,  it  is  still 
prohibitively  expensive.  To  bring  LES  closer  to  becoming  a  design  tool,  a  hybrid  RANS-LES 
approach  based  on  limited  numerical  scales  (LNS)  has  been  recently  developed  by  Metacomp 
Technologies  [12].  This  approach  combines  the  best  features  of  RANS  and  LES  in  a  single 
modeling  framework.  The  LNS  model  is  formulated  from  an  algebraic  or  differential  Reynolds 
stress  model,  in  which  the  sub-grid  stresses  are  limited  by  the  numerically  computed  local  length 
scale  and  velocity  scale  products.  LNS  thus  behaves  like  its  parent  RANS  model  on  RANS-type 
grids  but  reverts  to  an  anisotropic  LES  sub-grid  model  as  the  mesh  is  refined  locally,  thereby 
reaching  the  correct  (DNS)  fine  grid  limit.  Locally  embedded  regions  of  LES  may  be  achieved 
automatically  through  local  grid  refinement,  while  the  superior  near-wall  stress  predictions  of  the 
RANS  model  are  preserved,  removing  the  need  of  ad  hoc  topography-parameter-based  wall 
damping. 
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The  LNS  hybrid  formulation  is  well  suited  to  the  simulation  of  unsteady  flows,  including  mixing 
flows,  and  contains  no  additional  empirical  constants  beyond  those  appearing  in  the  original 
RANS  and  LES  sub-grid  models.  With  this  method,  a  regular  RANS-type  grid  is  used  except  in 
isolated  flow  regions  where  denser,  LES-type  mesh  is  used  to  resolve  critical  unsteady  flow 
features.  The  hybrid  model  transitions  smoothly  between  an  LES  calculation  and  a  cubic  k-e 
model,  depending  on  grid  fineness.  A  somewhat  finer  grid  was  placed  around  the  body  and  near 
the  jet,  and  the  rest  of  the  flow  field  is  occupied  by  a  coarser,  RANS-like  mesh. 

To  date,  the  LNS  technique  has  been  used  successfully  on  a  number  of  unsteady  flows. 

Examples  include  flows  over  cavities,  flows  around  blunt  bodies,  flows  around  airfoils  and 
wings  at  high  angle  of  attack,  separation  suppression  via  synthetic  jets,  forced  and  natural 
convection  flows  in  a  room,  and  mixing  flows  in  nozzles. 


3.  Results 


Time-accurate  unsteady  numerical  computations  that  use  viscous  Navier-Stokes  methods  were 
performed  to  predict  the  flow  field  and  aerodynamic  coefficients  on  the  subsonic  projectile  for 
“jet-on”  conditions.  Three-dimensional  numerical  computations  have  been  performed  for  the 
projectile  configuration  with  jet  interaction  using  CFD++  code  at  a  subsonic  Mach  number  of 
0.1 1  and  several  angles  of  attack  from  0°  to  4°.  The  preconditioned  version  of  the  CFD++  code 
was  used  to  obtain  efficient  numerical  solution  at  low  speed.  A  point-wise  two-equation 
turbulence  model  was  used  and  integrated  all  the  way  to  the  wall.  For  modeling  of  the  unsteady 
synthetic  jets,  a  hybrid  RANS-LES  approach  [12]  was  used.  For  computations  of  these  unsteady 
jets,  full  3-D  computations  are  performed  and  no  symmetry  was  used.  Computed  results  are 
compared  with  the  available  experimental  wind  tunnel  data  provided  GTRI2. 

For  computational  validation  purposes,  CFD  was  first  used  to  compute  the  flow  for  an  isolated 
2-D  jet  case  shown  in  Figure  1.  This  figure  shows  a  schematic  diagram  of  a  2-D  synthetic  jet 
along  with  a  flow  picture  obtained  from  the  experiment  [4].  Such  synthetic  jets  are  active 
control  devices  with  zero  net  mass  flux  and  are  intended  to  produce  the  desired  control  of  the 
flow  field  through  momentum  effects.  Here,  the  jet  width  is  0.5  mm  and  the  peak  jet  velocity  is 
20  m/s.  The  jet  actuator  operates  at  a  frequency  of  1000  Hz.  In  the  numerical  computations, 
unsteady  jet  boundary  conditions  were  applied  at  the  jet  exit  and  the  actual  cavity  was  not 
modeled.  A  sinusoidal  variation  was  used  for  the  jet  exit  velocity  with  a  peak  amplitude  of 
20  m/s.  Computed  velocity  and  vorticity  contours  are  shown  in  Figure  2.  The  time-averaged  jet 
centerline  velocity  over  many  cycles  is  compared  with  available  experimental  data  in  Figure  3 
and  is  found  to  be  in  reasonable  agreement.  As  shown  in  the  experiment  and  the  computations, 


2Private  communication  with  J.  McMichael,  Georgia  Tech  Research  Institute. 
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Centerline  axial  velocity  (m/s) 


the  time-averaged  centerline  velocity  is  decreased  with  increasing  distance  away  from  the  jet  exit 
(increasing  z  in  the  x-axis  in  Figure  3).  The  CFD  technique  described  earlier  was  then  applied  to 
compute  the  flow  over  a  projectile  at  subsonic  speeds. 


Figure  1.  Schematic  of  a  2-D  unsteady  jet  in  the  experiment. 


Figure  2.  Contours  of  velocity  components  (u,  v)  and 
vorticity  (from  left  to  right). 


Figure  3 .  Variation  of  time-averaged  centerline  jet  velocity 
with  distance  from  the  wall. 
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The  subsonic  projectile  is  a  1.8-caliber  ogive-cylinder  configuration  (see  Figure  4).  Here,  the 
primary  interest  is  in  the  development  and  application  of  CFD  techniques  for  accurate  simulation 
of  projectile  flow  field  in  the  presence  of  unsteady  jets.  The  first  step  here  was  to  obtain  the 
steady  state  results  for  the  same  projectile  without  the  jet.  The  converged  “jet-off”  steady  state 
solution  was  then  used  as  the  starting  condition  for  the  computation  of  the  time-accurate 
unsteady  flow  field  for  the  projectile  with  synthetic  jets.  Computations  were  also  performed  for 
the  steady  jet  cases.  The  jet  locations  on  the  projectile  are  shown  in  Figure  5.  The  jet  conditions 
were  specified  at  the  exit  of  the  jet  for  steady  (fixed  jet  velocity)  and  unsteady  (sinusoidal 
variation  in  jet  velocity)  jets.  The  jet  conditions  specified  include  the  jet  pressure,  density,  and 
velocity  components.  The  flow  field  inside  the  tiny  jet  cavity  is  not  computed.  For  the  unsteady 
jet  case,  time-dependent  jet  boundary  conditions  are  applied  at  the  jet  exit.  Numerical 
computations  have  been  made  for  these  jet  cases  at  a  Mach  number,  M  =  0.1 1,  and  at  an  angle  of 
attack,  a  =  0°.  The  jet  width  was  0.32  mm,  the  jet  slot  half-angle  was  18°,  and  the  peak  jet 
velocity  used  was  3 1  m/s  operating  at  a  frequency  of  1000  Hz. 


Figure  4.  Projectile  model  geometry. 


A  computational  grid  expanded  near  the  vicinity  of  the  projectile  is  shown  in  Figure  6.  Grid 
points  are  clustered  near  the  jet  as  well  as  the  boundary  layer  regions  to  capture  the  flow  regions’ 
high  gradients.  The  computational  grid  has  211  points  in  the  streamwise  direction,  241  in  the 
circumferential  direction,  and  80  in  the  normal  direction.  The  unsteady  simulation  took 
thousands  of  hours  of  central  processing  unit  time  on  a  Silicon  Graphics  Origin  computer 
running  with  24  processors. 
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Figure  6.  Computational  grid  near  the  projectile. 

Computed  pressure  contours  for  the  jet-off  ease  at  M  =  0.1 1  and  a  =  0°  are  shown  in  Figure  7. 
This  figure  shows  a  high  pressure  region  near  the  nose  and  low  pressures  in  the  near  wake. 
Computed  velocity  vectors  for  this  case  are  shown  in  Figure  8.  The  re-circulatory  flow  regions 
in  the  wake  and  near  the  step  upstream  from  the  base  are  clearly  evident.  As  seen  in  this  figure, 
the  flow  field  in  the  base  region  is  fairly  symmetrical.  The  flow  field  results  for  the  synthetic  jet- 
on  case  are  shown  in  Figures  9  and  10.  These  figures  show  the  qualitative  flow  features  near  the 
jet  as  well  as  the  base  region  of  the  projectile.  Figure  9  shows  the  velocity  vectors  at  a  given 
instant  in  time.  It  clearly  shows  the  flow  in  the  base  region  to  be  asymmetric  because  of  the 
interaction  of  the  unsteady  jet.  The  shear  layer  from  the  free  stream  flow  resulting  from  the  step 
corner  upstream  from  the  base  interacts  with  the  unsteady  jet  and  deteriorates  just  a  short 
distance  downstream  from  the  jet.  The  jet  has  a  strong  effect  on  the  pressure  distribution 
upstream  and  downstream  from  the  jet  location  (see  Figure  10).  The  pressure  field  upstream  and 
downstream  from  the  jet  is  clearly  affected  by  the  jet  flow,  depending  on  whether  the  flow  is 
entering  or  leaving  the  cavity.  The  computed  flow  field  again  is  asymmetrical.  Figure  10  also 
shows  regions  of  alternating  low  and  high  pressure  just  downstream  from  the  jet,  indicating 
vortieity  being  shed  from  the  jet  exit  into  the  base  region  flow. 

The  resulting  surface  pressures  from  the  unsteady  flow  fields  are  integrated  to  obtain  the 
aerodynamic  forces  and  moments.  The  computed  axial  force  and  normal  force  are  shown  in 
Figures  1 1  and  12  as  a  function  of  time,  respectively.  These  computed  results  clearly  indicate 
the  unsteady  nature  of  the  flow  field.  These  forces  are  found  to  change  as  a  function  of  time. 
Changes  in  the  drag  force  are  smaller  than  those  of  the  lift  force  over  the  first  50  ms.  In  both 
figures,  computed  results  obtained  with  unsteady  RANS  (URANS)  as  well  as  the  hybrid  RANS- 
LES  approach  referred  to  here  as  the  LNS,  are  compared  to  each  other.  As  shown  in  Figure  1 1, 
LNS  and  URANS  both  produce  a  non- zero  lift  force  attributable  to  the  jet  at  0  degree  angle  of 
attack.  The  lift  computed  by  the  URANS  approach  is  about  half  that  computed  by  the  LNS 


approach.  As  shown  in  Figure  12,  the  drag  force  variation  as  a  function  of  time  is  small  with  the 
LNS  method,  whereas  the  variation  is  almost  negligible  with  URANS.  Although  not  shown 
here,  the  pitching  moment  also  shows  a  variation  similar  to  that  of  the  lift  force. 


Figure  7.  Computed  pressures,  jet  off,  M  =  0.11,  a  =  0°. 
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Figure  1 1 .  Computed  lift  force,  M  =  0. 1 1 ,  a  =  0°. 
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Figure  12.  Computed  drag  force,  M  =  0.1 1,  a  =  0°. 

Figures  13  and  14  show  the  time  histories  of  the  lift  and  drag  forces,  respectively,  for  various 
angles  of  attack.  These  computed  results  are  obtained  with  the  LNS  turbulence  model. 

Figure  13  shows  the  increase  in  the  lift  force  resulting  from  the  jet  interaction  for  all  angles  of 
attack  considered.  The  values  for  the  jet-off  condition  correspond  to  time  zero.  As  shown  in 
Figure  14,  the  drag  force  is  also  seen  to  increase  somewhat  because  of  the  jet  interaction,  but  the 
effect  of  angle  of  attack  on  the  drag  force  is  rather  small. 

The  unsteady  interaction  of  the  synthetic  jet  and  the  wake  flow  of  the  projectile  usually  require 
hundreds  of  cycles  of  the  jet  operations  before  a  truly  unsteady  periodic  flow  is  established. 
Computations  were  extended  to  50  additional  cycles  of  jet  operation  for  0-  and  4-degree  angles 
of  attack.  These  computed  results  are  shown  in  Figure  15  as  a  function  of  time  until  100  ms. 
Also,  in  order  to  test  the  time  accuracy,  the  time  step  was  reduced  from  0.025  ms  to  0.00833  ms 
approximately  at  time  70  ms.  Although  not  shown  here,  the  number  of  inner  time  steps  in  the 
dual  time-stepping  scheme  used  here  was  also  varied  and  the  results  were  found  to  remain 
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essentially  unaltered.  As  can  be  seen,  the  mean  value  of  the  lift  force  for  0-  and  4-degree  angles 
of  attack  did  not  change  much.  For  the  0-degree  angle  of  attack  case,  the  mean  value  of  the  lift 
force  was  calculated  to  be  about  0.18  Newton,  resulting  from  the  unsteady  jet  interaction. 
Similarly,  time-averaged  mean  values  of  the  lift  forces  at  other  angles  of  attack  (a  =  2  and  4) 
were  computed.  The  difference  between  these  computed  mean  values  of  the  lift  force  and  the 
corresponding  lift  at  jet-off  condition  then  yields  the  net  change  in  the  lift  force  attributable  to 
the  jet  interaction.  Figure  16  shows  the  change  in  lift  force  attributable  to  the  jet  interaction  at 
various  angles  of  attack.  The  computed  results  here  are  found  to  be  generally  in  good  agreement 
with  the  experimentally  [17]3  obtained  lift  force  attributable  to  the  jet.  As  seen  in  both  CFD  and 
the  experiment,  the  change  in  lift  attributable  to  the  jet  interaction  is  almost  constant  with  angle 
of  attack  from  0°  to  4°. 

1.2  | - - - 1 

1  ' 

08  i 

0.6  A 

Fy  0.4 

0.2 

0 

.0  2  -  - LNS-A2 

- LNS-A4 

-0.4  ' - 1 - ' - 1 - 1 _ 

0  10  20  30  40  50 

Time  (ms) 

Figure  1 3 .  Computed  lift  force  for  various  angles  of  attack,  M  =  0. 1 1 . 
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Also,  a  large  number  of  CFD  data  sets  were  saved  at  regular  intervals  to  analyze  the  fully 
unsteady  periodic  nature  of  the  flow  field  which  required  huge  computer  resources  in  terms  of 
storage  and  flow  visualization.  The  computed  results  obtained  thus  far  have  been  used  to 
determine  the  feasibility  of  the  synthetic  jets  to  provide  control  authority.  The  next  challenge  is 
to  use  the  current  methods  with  the  same  projectile  but  with  spin.  Future  efforts  will  require 
large  unsteady  numerical  computations  to  optimize  the  number  and  locations  of  the  synthetic 
jets.  It  is  anticipated  that  multiple  jets  may  be  required  to  provide  the  control  authority  needed 
for  maneuvering  a  subsonic  projectile  and  will  require  the  use  of  a  coupled  CFD  and  controls 
technique  to  divert  control  authority  by  means  of  a  jet. 


3Private  communication  with  J.  McMichael,  Georgia  Tech  Research  Institute. 
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Figure  14.  Computed  drag  force  for  various  angles  of  attack,  M  =  0.1 1. 
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Figure  15.  Computed  lift  force  for  various  angles  of  attack, 

M  =  0.1 1,  for  a  longer  time  period. 


Figure  16.  Computed  change  in  lift  force  attributable  to  jet  at 
various  angles  of  attack,  M  =  0.1 1,  a  =  0°. 


4.  Concluding  Remarks 


This  report  described  a  computational  study  undertaken  to  determine  the  aerodynamic  effect  of 
tiny  synthetic  jets  as  a  means  to  provide  the  control  authority  needed  to  maneuver  a  projectile  at 
low  subsonic  speeds.  Computed  results  have  been  obtained  at  a  low  subsonic  speed  and  0  degree 
angle  of  attack  for  a  subsonic  projectile  via  the  Navier-Stokes  computational  technique  and 
advanced  turbulence  models.  Qualitative  flow  field  features  show  the  interaction  of  synthetic  jet 
with  the  free  stream  flow  and  the  large  extent  of  influence  upstream  and  downstream  from  the  jet. 
Unsteady  jet  results  obtained  for  a  2-D  jet  are  compared  with  the  data  and  are  found  to  be  in 
reasonable  agreement.  The  unsteady  jet  in  the  case  of  the  subsonic  projectile  is  shown  to  sub¬ 
stantially  alter  the  flow  field  near  the  jet  and  the  base  region,  which  in  turn  affects  the  forces  and 
moments  even  at  0  degree  angle  of  attack.  RANS  and  hybrid  RANS-LES  turbulence  models  were 
used  in  the  computations  at  0  degree  angle  of  attack.  The  lift  predicted  by  the  RANS  simulation  is 
less  than  half  the  predicted  lift  with  the  hybrid  model.  The  predicted  changes  in  lift  force  attribu¬ 
table  to  the  jet  match  well  with  the  experimental  data  for  various  angles  of  attack  from  0°  to  4°. 

The  results  show  the  potential  of  CFD  to  provide  insight  into  the  jet  interaction  flow  fields  and 
provide  guidance  as  to  the  locations  and  sizes  of  the  jets  to  generate  the  maximum  control 
authority  for  maneuvering  smart  munitions.  Current  and  future  efforts  are  directed  toward  the 
addition  of  spin  and  yaw  to  these  large  unsteady  synthetic  jet  computations  for  application  to  a 
spinning  projectile.  Future  efforts  will  also  include  the  optimization  of  the  number  and  location 
of  the  synthetic  jets  as  well  as  coupling  of  CFD  and  controls  to  accurately  determine  the  control 
authority  needed  for  maneuvering  a  subsonic  spinning  projectile. 
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